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O ■ ABSTRACT 

Diskoseismology, the theoretical study of small adiabatic hydrodynamical 

global perturbations of geometrically thin, optically thick accretion disks around 

■ black holes (and other compact objects), is a potentially powerful probe of the 

. gravitational field. For instance, the frequencies of the normal mode oscillations 

On ' 

can be used to determine the elusive angular momentum parameter of the black 
t:J- ■ hole. The general formalism developed by diskoseismologists for relativistic sys- 



tems can be readily applied to the Newtonian case of cataclysmic variables (CVs). 
O ■ Some of these systems (e.g., the dwarf nova SS Cygni) show rapid oscillations in 

the UV with periods of tens of seconds and high coherence. In this paper, we 
^ ■ assess the possibility that these dwarf nova oscillations (DNOs) are diskoseismic 

■ modes. Besides its importance in investigating the physical origin of DNOs, the 

present work could help us to answer the following question. To what extent 
are the similarities in the oscillation phenomenology of CVs and X-ray binaries 
(XRBs) indicative of a common physical mechanism? 

Subject headings: accretion, accretion disks — hydrodynamics — stars: indi- 
vidual (SS Cygni) — stars: individual (VW Hyi) — stars: novae, cataclysmic 
variables — white dwarfs 
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Introduction 



During the outburst pliase (timescale ~ a few days) some non-magnetic {B <^ 10^ G) 
CVs exhibit fluctuations with periods P ^ 10 s and Q = l/\dP/dt\ ~ 10^ - 10^. The 
phenomenology of these DNOs is very rich (jWarnerll2004l ). Also observed are IpDNOs (long 
period DNOs): P > 30 s and Q ~ 10^ — 10^, as well as QPOs (quasi-periodic oscillations): 
P > 100 s and Q ~ 10. At times, all three features can be observed simultaneously. 

The fact that P ~ the Keplerian frequency at the radius of the white dwarf has moti- 
vated several candidate theoreti cal explanations for DNOs, ranging from 'hot blobs' (jPopham 
19991 ) to star oscillation modes ( jPapaloizou &: Pringlelll978l ) to explanations based on the dy- 
namics of the boundary layer or some other type of transition region located where the disk 
touches the white dw 'arf . Examples of the las t type include magnetic accretio n onto a slipping 
belt at the equator (IWarner fc Woudtll2002l ) and 'spreading layer' models (jPiro &: Bildsten 
2004r). Hvdrodynan iic os cillations in and near the boundary layer have been investigated by 
Carroll et all Jl985h and ICoUins. Heifer fc Van koO(k_ However, their analyses were 

local (plane-wave WKB) in radius and height. Indeed, ICarroU et al. ( 1985 ) mention the 
need of a global analysis. None of these models is truly successful at the quantitative level, 
although magnetic models appear to be favored. 

Lately, o ther proposals based o n the dynamics of the accretion disk have appeared. 
For example, iKluzniak et al.l (120051 ) propose a nonlinear hydrodynamical disk resonance 
explanation for the oscillations. Part of the reason this approach is attractive is that it could 
explain the similarity of the phenomenology in CVs and XRBs. 

Some DNOs exhibit a 1:2 (e.g., the period halving in SS Cyg) and in one case (VW Hyi) 
a 1:2:3 harmonic stru cture. In sorne sources, a 1:15 relatio n of the frequencies of a QPO 
and DNO is observed jMauchel I2OO2I : IWarner fc Woudtl l2005h , similar to that seen in XRBs. 
Some high-frequency QPOs in black hole XRBs are in a 3:2 ratio (lAbramowicz fc Kluzniak 
2OOII : iMcClintock fc Remillardll2006l ). All the explanations listed above (except the hot blob 
and resonance models) are not applicable to black holes, given the absence of a surface (even 
though some could apply to neutron stars). Conversely, QPO explanations based on general 
relativistic effects are not applicable to the CV case. The only common explanation is based 
on hydrodynamical oscillations in the accretion disk (and/or corona). 

With the intention of adding to this debate, this paper explores the possibility that the 
DNOs in CVs are caused by global hydrodynamical oscillations (normal modes) in the inner 
accretion disk. To this end we will use diskoseismology, a formalism previously developed 
to stu d y small adiabatic perturbations in relativ i stic a ccretion disks [see, e.g., [Wagoner 
( I1999I ): IWagoner. Silbergleit. fc Ortega- Rodriguea (I2OOII )]. The DNOs would arise as the 
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fluctuations modulate the outgoing radiation. Although lYamasaki. Kato &: Mineshigd (Il995l ) 
mention global hydrodynamical p-modes (trapped acoustic waves) in accretion disks around 
white dwarfs, these modes lie in a different region from the one studied in this paper; in 
addition, the physical analysis of these authors is much simpler than ours, as they work with 
vertically integrated variables and limit themselves to axisymmetric modes. 

In section 2 we briefly discuss the unperturbed disk model and then outline diskoseis- 
mology in the general case. The rest of the paper apphes diskoseismology to cataclysmic 
variables. An introduction (section 3) precedes the formal WKB solution in section 4. Fi- 
nally, section 5 presents the results and discusses their implications. 



2. Diskoseismology 

We take c = 1, and express all distances in units of GM/(? and all frequencies in units 
of (? jGM (where M is the mass of the white dwarf) unless otherwise indicated. We work 
in cylindrical coordinates r,ip,z. 

The stationary [d/dt = 0), symmetric about the midplane z = 0, and axially symmetric 
{d/dip = 0) unp erturbed disk is taken to b e described by the sta ndard thin disk models for 



the Newtonian (IShakura fc Sunyaevi Il973l ) and the relativistic (iNovikov &: Thornd Il973l ) 
cases. In particular, we note that the velocity components = = (the disk fluid 
moves in nearly circular orbits), the disk semi-thickness h{r) ~ Cs/f2 <^ r (where Cs and 
Q are the speed of sound and the Keplerian frequency, respectively), the disk is optically 
thick and radiatively efficient, the viscous stress is given by ap (i.e., the pressure multiplied 
by a constant parameter), and typical values for the disk density and central object mass 
make the neglect of self-gravity, at both the unperturbed and the perturbed level, a good 
approximation. Further properties of the unperturbed disk will be discussed in section 5. 

Akin to helioseismology, diskoseismology is the theoretical study of the normal modes 
of small adiabatic hydrodjTiamical global perturbations in geometrically thin, optically thick 
accretion disks. In this section, we describe diskoseismology in general terms. In the following 
sections, we will use Newtonian diskoseismology to study oscillations in white dwarf accretion 
disks. Fortunately, the relativistic and nonrelativistic models are formally very similar. 

Even though the unperturbed disk is not in equilibrium, this perturbative approach is 
justified whenever the viscous time scales during which the disk evolves as a whole are much 
longer than the modes' periods, and whenever the turbulence has spatial scales which are 
smaller than the modes' spatial scales, and it is quasi-steady on the time scales of the modes. 
(These conditions hold for the modes dealt with in this paper.) 
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We consider barotropic disks [p = p{p), giving vanishing buoyancy frequency] and as- 
sume that the results thus obtained are generic, since typically the buoyancy frequency is 
much less than the characteristic dynamical frequency Q. In this case hydrostatic equilibrium 
provides the vertical density and pressure profiles 



p = Po{r){l - y y , p = PQ{r){l-y) 



2\g+l 



g = l/iT-l)>0 



where F > 1 is the adiabatic index. One has F = 4/3 within any radiation pressure dominated 
region of the disk, and F = 5/3 within any gas pressure dominated region (such as the white 
dwarf accretion disks). The disk surfaces are at y = ±1, with y related to the vertical 
coordinate z by 

z 

y 



h{r) 



1 



2F 



To investigate the eigenmo des of the disk osc i llations star ting from these assumptions, 
one applies the formalism that llpser &: Lindblonj (jl99ll . Il992l ) developed for perturbations 
of purely rotating perfect fluids, starting from the equations of conservation of mass and 
momentum. One can thus express the Eulerian perturbations of all physical quantities 
through a single function 6V oc Sp/p (the Eulerian perturbation of the pressure divided by 
the density) which satisfies a second-order partial differential equation. Due to the stationary 
and axisymmetric background, the angular and time dependences are factored out as SV = 
V{r, z) exp[i{rrup + crt)], where a is the eigenfrequency, and the governing PDE becomes 



d_ 

dr 



9"-p 



dV 



— dr 



+ 



d_ 

dz 



p dV 

uS^ dz 







(2) 



where /3 = dt/dr, g^^ is the metric coefficient, Cg stands for the speed of sound, the corotation 
frequency 

uj = a + mQ , (3) 

Q{r) is the Keplerian frequency, and K,{r) is the radial epicyclic frequency (i.e., the frequency 
of radial perturbations of free particle circular orbits). 

An assumption of strong variation of modes in the radial direction (characteristic radial 
wavelength <^ r) ensures the approximate WKB separability of variables in this PDE, 
V{r,z) = Vr{r)Vy{r,y). T he 'vertical' part , Vy, of the functional amplitude V{r,z) varies 
slowly with r. They obey (jPerez et al.l 119971 . hereafter referred to as RDl) 



d'^Vr 

dr'^ 



[i-y 
1 



dy^ 
' d , 



dV,. n 

2gy-^ + 2gul 



dy 



1- 1 



(CU^ — K^) 



dr 



2 2\ 



dVr 

dr 



K = 



Vr = 



(4) 
(5) 
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where 



= uj/Q± ; 



(6) 



Qx(r) refers to the vertical epicychc frequency (i.e., the frequency of vertical perturbations 
of free particle circular orbits), and the speed of sound is evaluated at midplane. In these 
equations ^(r) is the slowly- varying separation function. 

These equations have been used in the past to reveal the properties of different types of 
modes in relativistic accretion disks. In this paper we solve the equations in the Newtonian 
limit. We employ, as in the original papers, WKB methods in a straightforward manner. 

3. Diskoseismic Modes for the White Dwarf Case: Physical Parameters and 



We now apply diskoseismology to the accretion disks of white dwarfs. We are thus 
working in the Newtonian hmit (for a mass M — Mq white dwarf, its radius i? ~ 5 x 10^, 
corresponding to 7.4 x 10^ km), which amounts to setting g^'^ — P — 1, g — 3/2 in the 
equations of the previous section. We note in particular that 



i.e., the rotational, radial epicyclic, and vertical epicyclic frequencies are all equal (up to 
corrections of relative size ~ 1/r or smaller). As we will see, this results in a qualitative 
change in the nature of the solutions with respect to the relativistic case. 

This is a straightforward problem with few free parameters. In addition to the white 
dwarf mass Af, there is [> R{M)], the inner radius of the accretion disk. Then, one 
needs a functional form for the speed of sound. It is fortunate that the properties of the 
unperturbed disk (viscosity and luminosity) enter into the formalism only through the speed 
of sound, and that the dependence of the speed of sound on these properties is rather weak. 

In this exploratory calculation, we ignore the complications related to the nature of the 
transition from the disk to the star surface or magnetosphere and assume that DNO oscil- 
lations are confined to the disk. We introduce a disk inner boundary condition parameter, 
given our ignorance of the physical conditions in the (magnetic) boundary layer. As we will 
see, the results are rather insensitive to this parameter. The approximate character of our 
analysis should be stressed. A deeper study of boundary layer dynamics is needed in order 
to understand better the physical conditions that would allow for the existence of the modes 
described in this paper. For instance, one could ask how significant is the damping and 
'leakage' of the modes into the magnetosphere, where the radial epicyclic frequency is given 



Scope 




(7) 
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by the more general expression = 2Qr ^d{r'^Q)/dr, if the modes are still confined to a 
small vertical extent h <^r. 

As in the case of relativistic diskoseismology (with black holes or weakly magnetic 
neutron stars), the formalism reveals the existence of different types of modes. In this first 
diskoseismic investigation of cataclysmic variables, we choose to study the g-modes, since 
they are trapped near the inner edge of the disk and are therefore the most luminous (and 
robust). By definition, the g-modes have corotating eigenfrequencies which are smaller than 
the gravitational radial oscillation (epicyclic) frequencies at the involved radii. In other 
words, the (smaller) pressure restoring forces act against the gravitational restoring forces. 
(The 'g' in the name comes from a mathematical similarity to the g-modes within slowly 
rotating stars.) The acoustic (p-) modes extend to the outer region of the disk, thus involving 
strong effects of viscosity and the uncertain boundary conditions at its outer radius. 



4. Formal WKB Solution of the Eigenvalue Problem for G— modes 



In this section we will proceed to solve formally the eigenvalue problem in order to 
obtain equation (|28|1 . our 'master' equation, from which all further results in the paper are 
derived. 



Taking into account expressions ([7j), the separated equations (jlj) and 
Newtonian limit to 



reduce in the 



d^Vr 

dr"^ 



1 



d'^Vy 

dy"^ 



Q?) 
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O dVy 
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i^-y') 









dVr 

dr 



K = 



K = 



(9) 



(Recall from section 2 that (7 = 3/2 for white dwarf accretion disks.) Equations ([8]) and ([9]) 
are two coupled eigenvalue equations, the eigenvalues being oj and \E' (or equivalently, a and 
'^) . We will solve these equations within a straigh tforward WKB formalism as described 
by lOrtega-Rodrfguez. Silbergleit. fc Wagonerl (120021 ) (hereafter referred to as RD3). With 
respect to boundary conditions, it turns out that it is sufficient to require that the vertical 
eigenfunction and its derivative remain finite at the surface of the disk and to introduce an 
inner boundary parameter for the radial equation, as mentioned. 

By definition, g-modes exist wherever uj"^ < /t^(= Vt^). This means that the mode exists 
between an inner radius (which is r^, the inner radius of the accretion disk, for the cases 
studied in this paper) and the radius r+, defined as the radius at which \uj\ = Q. The mode 
decays exponentially for r > r+. 



- 7- 



4.1. WKB Solution of the Radial Eigenvalue Problem 



We will cast equation ([9]) in a form amenable to a WKB solution, using a procedure 
which is similar to the one developed in RD3. We introduce the dependent variable 

' (10) 



W 



— LU^ dr 



and the independent variable 



T 



dr' 



so that equation ([9]) can be written in the simpler form 



(12) 



The mode exists in the range ri < r < r^ or, equivalently, < r < = 'r(rj). The function 
S{t) changes sign at r = (where r = r+). The boundary conditions are that W decays 
for r < 0, while at the other extreme, at r = r^, we parameterize our ignorance at the inner 
boundary by means of the quantity 9i: 



(n"^ -co^) cos 9iW -sin 9, 



dW 
dr 







which is equivalent to 



cos 9i—--^ — sin 9iVr = 
dr 



(13) 
(14) 



A straightforward WKB solution can be obtained for equation f|T2l) at all points except for 
r close to 0: 



W oc S~^/^(t) cos 



(15) 



where $ is a constant to be determined. For r near zero, the solution is given by the Airy 
function of the first kind: 



W oc Ai{-s]l'^T) , s+ = dS{0)/dT > . 



(16) 



The asymptotic matching of solutions f[T^ and f[T^ in the usual manner yields $ = Tr/l+nvr, 
where n is an integer, and 



c,\r)^/{n^ - cu2)(^/cu2 - l)dr = Tr{n ± 1/4) 



(17) 



for the cases 9i = 7t/2 and 0, respectively. Thus n is the radial mode number constrained by 
the fact that the integral in equation f|T7j) is positive. 
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4.2. WKB Solution of the Vertical Eigenvalue Problem 

Using the new independent variable 



we can write equation ([8]) in the simpler form 



which is amenable to a WKB solution. 



(19) 



We now discuss the issue of boundary conditions. In view of the symmetry of the 
differential equation ([8]), one has 



dy 







or 



y=0 







(20) 



for even and odd modes in the disk, respectively, and it is enough to consider the interval 
< y < 1. Near the singular boundary point y = I, equation (IHl) can be written 



3 dy^ 



1 



+ 



dyy_ 

dy 



+ 



+ . . . K = 



(21) 



where the dots repres ent no n-singular terms. According to the analytical theory of second 
order ODEs [see, e.g., Olver (fl982 )]. the general solution to equation ([8]) near y = 1 has the 
form 

Vy = CMy) + C2(i - yV^vM , (22) 

where Ci and C2 are arbitrary constants, and Vi{y) and V2{y) are converging power series of 
(1 — y). Imposing the boundary condition that the eigenfunction and its derivative be finite 
aX y = 1 implies thus that C2 = 0. 

As noted in RDl and implied by equation ( 12T1) and the boundary condition, the vertical 
eigenfunction near the boundary can be written in terms of the Bessel function of the first 
kind as 

Vy{y) cc (l-y)-^/Vi/2(2K|v/(3/2)(l-y)) , (23) 
the asymptotic (large argument) expression of which is 



Vyiy) oc (1 - y)-i/2cos [|^,|v/6(l-y)- (7r/2) 



(24) 
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On the other hand, equation (fT9l) can be solved with a WKB method for values of y not 
close to 1: ^ 

ocQ-^/^(r)cos I Q^'^{C')dC' -1tx/2 , (25) 



where T = 0, 1 for the even and odd eigenfunctions, respectively. Matching this with the 
asymptotic expression of the boundary solution (12^ . we arrive at 




l + (vl//a;2-l)(l-y2) 



dy 



ttJ 



id* 



(26) 



The parity (of the vertical eigenfunction) and the vertical mode number j are contained in 
a single quantity J = j + 3/4±l/4for the odd and even cases, respectively. Thus J takes 
positive integer and half- integer values (see Table 1). 

It should be stressed that, because of the method used to solve the equations, the above 
results are accurate only when "^/wl — 1^1 for all r in the range < r < r+. 



4.3. Uncoupling 

Equations (fTTl) and (l26l) are two eigenvalue equations (the eigenvalues being \1^ and u), 
which need to be uncoupled. We start by inverting the integral / on the left hand side of 
equation (|26l) to find \E'/u;^ — 1 as a function of /: 

^/col - 1 = /(/) . (27) 

Once this is done numerically, the function /(/) can be substituted in equation ( ITTl) : 

/ ^ c;l(r)■^(^]2 _ uj^)f{nJ/V3\u,\)dr = n{n ± 1/4) . (28) 

The left-hand side of equation (1251) now only depends on the eigenfrequency a and not on 
\l/ [recall definition ([3])]. All further results in this paper are obtained from this equation. 



5. Results and Discussion 

In order to use equation fl28l) and look for solutions with integer values of n > 0, we 
need a value for and a functional form for the speed of sound Cg. We will take = 5000, 
corresponding to 73S3{M / Mq) km. The restriction rj > R, the white dwarf radius, places a 
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limit on its mass M. Using a standard mass-radius relation (jNauenberg||l972l ). this value of 
Ti implies that M > 0.76Mq. 



As for the speed of sound, there are different solutions given by IShakura fc Sunyaev 



(I1973I ). depending on the region of the disk in which one is working. For the values of 
interest in this paper, the relevant expression is the one corresponding to the outermost 
solution, in which temperatures are low enough that pressure is dominated by gas (not 
radiation) and opacity is dominated by photon absorption in the form of nonrelativistic 
thermal bremsstrahlung ('free-free' transitions). 

Their solution for the speed of sound is given by 

(cjc)-^ = 300 a^/^°(M/M0)^/^M-'/'° r^^^ , (29) 

where c, a, M, and Miq stand, respectively, for the speed of light, the (dimensionless) 
viscosity parameter, the mass of the white dwarf, and the mass accretion rate in units of 
10^^ g/s. The small powers of the three parameters render our results relatively insensitive 
to them. We remind the reader that the interior structure of the unperturbed accretion 
disk enters our formulation only through the speed of sound. As mentioned above, we 
parameterize our ignorance about the radial inner boundary by means of the quantity 6i 
introduced in equation (fn|) . 



5.1. Axisymmetric Modes 



The eigenmodes are described by the azimuthal, radial and vertical mode numbers m, 
n, j, plus the value of the parity of the vertical eigenfunction. We start by studying the 



m = modes. For them, the condition Q"^ —u 
the modes are trapped between and r_|_ = |cr 



^ > becomes fl^ 

2/3 



cr^ > 0, which means that 



Using the characteristic values = 5000, M = Mq, a = 0.1 (lLasotall200ll ). and Mig = 1 
(IFrank. King fc Raindl2002l ). we calculate a for different values of n, j and parity for the 
case 6i = 7r/2. (Note that the values of j and parity are contained in J in the manner 
described in section 4.2 and Table 1.) Table 2 shows the results. The dimensional frequency 
/ = 3.23 X 10^{M(7)/M)\a\ Hz. The n = 0, J = 1 entry thus corresponds to a period of 12 
seconds. 

Table 3 shows Ar/vi = (r_|_ — ri)/ri for the same values of the parameters. As in the 
case of axisymmetric g-modes in relativistic accretion disks (RDl), the value of c^/c is small 
enough that: a) the small n eigenfrequencies are close to the maximum value of n{r) [which 
equals f2(rj) = 2.83 x 10~^ in the present case], approaching it as j increases; b) the modes' 
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radial width decreases as j increases and increases as n increases. 

Self-consistency requires that the parameter "^/ool — 1 not be small, and the larger it 
is, the better the WKB approximation. It turns out that \l'/ci;^ — 1 ~ 1 when J = 1 and 
\E'/ct;^ — 1 ^ 1 when J > 1, so the procedure passes this test. (Incidentally, it is interesting 
that the minimum value of \E'/ct;^ — 1 over the allowed values of the radius does not depend 
on n or m, only on J.) 

There is no J = 1/2 entry in Table 2 because there is no WKB solution. A direct 
numerical solution of equations ([8]) and iQ shows that the absence of a solution is not 
an artifact of the WKB method, although it could result from the assumptions made in 
obtaining these equations. The possibility remains open that there are modes not detected 
by our methods. 

Next we compute the behavior of the eigenfrequency a as the various parameters change. 
For the choice J = 1, n = 0, 6i = 7r/2, we obtain 

d\oga/d\ogri = —1.5 , (30) 

9 log a/9 log a = 0.0041, (31) 

aioga/aiogM = 0.010, (32) 

aiogcj/aiogM = -0.0061, (33) 

which means that the dependence of a on a, M, and M can be neglected since it is much 
weaker than its (expected) dependence on rj. Moreover, a change of di from 7r/2 to implies 
only a change of a few percent in a. Note, however, that there can be no n = solution 
with 9i = because then the RHS of equation (fT71) would be negative. 

Also note that equation (1321) assumes a fixed value of r^. Care must be taken then 
when comparing with observations, since usually the parameters are not independent. For 
example, for a set of white dwarfs of fixed composition there is a relationship between their 
radius (R) and mass. So since a = cr(ri) and = ri{M,M) (for fixed surface magnetic 
field), one has thus, for fixed M, 

dlogM dlogM dlogri ' ' ^ ^ 

for the case M oc r^^ f lNauenberg Il972l ) appropriate to low-mass white dwarfs with weak 
magnetic fields (r^ = R). 

When computing the variation of a with the luminosity L for any specific cataclysmic 
variable, if one assumes that oc L~'^, where 7 is a constant, one obtains 

d log P d log r j d log a 



d log L d log L d log r j 



-1.57. (35) 
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Such a dependence of on luminosity could be due to strong magnetic field effects (IShapiro &: Teukolsky 



19831 ). The idea is that inside the magnetosphere the accretion flow is controlled by magnetic 
fields, so one takes rj to be the equatorial radius of the magnetosphere. Its value should be 
only moderately greater than the white dwarf radius R for those magnetic field strengths 
thought to be present in many of those white dwarf systems which exhibit DNOs. F or exam- 



ple, th e value of d\ogP/d\ogL = —0.091 for the dwarf nova SS Cygni reported by iMauche 



( 1l996l ) is reproduced by our formalism provided 7 = 0.06. 



In summary, the model may be consistent with observations. The radial width of the 
modes may be large enough to make them visible, and the periods are consistent, although 
this is hardly surprising as they are related to the Keplerian period at rj. This leads to 
the limit P > lO.94(M0/M)^/^[rj(km)/7383 km]^/^ seconds on the periods of axisymmetric 
modes. However, the physics of the magnetosphere is at present too uncertain to provide a 
means to test the model accurately. At this point it makes sense to look at the higher m 
modes. 



5.2. Non- Axisymmetric Modes 

Here we discuss the m > modes. Note that from definition ([3]) and the fact that only 
|a;| enters in the equations, we can, without loss of generality, consider only the non- negative 
m cases provided we allow for negative a. 

The condition — a;^ > implies that 

-m -1 < a/n{r) < -m + 1 . (36) 

Therefore, g-modes are trapped between and r_|_ = (m+ l)^/^|cr|^^/^. The inner boundary 
of the mode (r_) is r^, provided that \a\ > {m — l)Q{r) . The negative sign of a reflects the 
fact that the oscillations' wavefronts travel in the same direction as the fluid. 

Table 4 and Figure 1 show some properties of the n = 0, J = 1 modes for the m = 0, 1 
and 2 cases. Interestingly, the results are consistent with the obser vation of a 1:2:3 harmoni c 



structure in the power spectra of the DNO VW Hyi, reported by IWarner fc Woudtl (120051 ). 
The observed peaks remain harmonic even during a threefold increase in the fundamental 
frequency. As in other models, the increase in frequency is due to a decrease in induced 
by an increase in M (and thus L). 

If one expresses a in terms of r_|_ and works to lowest order in Ar/r^, from equation (l28l) 
one obtains 

Ar/vi ~ (n ± 1 /4)2/3(m + 1)-^/^^^ . (37) 
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Since e = h{ri)/ri ~ (c^/fir),.., using equation fl29|) and our chosen set of parameters (which 
only weakly affect the estimate) gives e ~ 4.2 x r\^^ ^ 1.2 x 10"^. Then equation fl37j) 
is seen to agree approximately with the values in Tables 3 and 4. 

The m > diskoseismic modes are in principle harder to observe than the axisymmetric 
modes, since in the former the modulation of the outgoing radiation is not as efficient due to 
cancelations over the entire disk (and their radial extent is somewhat smaller). However, the 
fact that the white dwarf can eclipse parts of the inner disk renders these modes more visible. 
We predict that the energy spectrum of the DNOs is characteristic of the very inner disk 
(and not other regions such as the magnetosphere), and thus should be close to blackbody 
at the highest disk temperature. 

As in the case of g-modes in relativistic accretion disks, the eigenfrequency splittings 
when m changes by unity are much larger than the n and j splittings. The g-modes that 
are trapped by the general relativistic behavior of the radial epicyclic frequency k have 
m splittings that are not as harmonic as those found here (Table 4). In particular, their 
eigenfrequencies are only proportional to m for m > 2 [an interesting difference between the 
relativistic and the nonrelativistic eigenfrequencies is that, for m ^ 1, |cr| mQ{ri) for the 
former (see RDl) while |cr| (m + l)Q{ri) for the latter]. However, it must be noted that 
it appears that g-modes with small mode numbers are viscously unstable, at least if the 
vis cosity acts in the usual hydrodynami c way. This result is proven for the relativistic case 



m 



Ortega-Rodriguez &: Wagoner! (120001 ). where the exponential growth rate of the mode 
is roughly afl. For the nonrelativistic g-modes studied in this paper, an application of 
their method yields the same result. Furthermore, these low-lying modes do not contain a 
corotation resonance [where u{r) = 0]. 



Comments 



Thus, if the ideas in this paper are correct and DNOs are linear diskoseismic modes 
then the results could suggest a scenario in which white dwarf DNOs and b lack hole QPOs 
are pr oduced through somewhat different physical mechanisms. After all, as iKluzniak et al. 
( I2OO5I ) themselves point out, the phenomenology is similar but not identical. While in 
white dwarf DNOs the simultaneous observation of two or even three of the harmonics is 
well established, in black holes usually only one mode of oscillation is observed at a given 
moment. Furthermore, white dwarf DNOs (and neutron star QPOs) have frequencies that 
vary with time, while QPOs in black holes have fixed frequencies. In neutron star XRBs, the 
(twin) high-fre quency QPOs are analogous to the DNOs, but their frequency ratio is much 
closer to unity (Ivan der Klisll2006l ). It appears to be controlled (magnetically) by the neutron 
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star, since their frequency difference is close to the neutron star spin frequency or half of it. 
We note that some DNOs are split at a QPO frequency. We also note that with increasing 
radius of a neutron star magnetosphere, their g-modes should approach (for > 8) those of 
this paper. 

We hope to investigate the magnetic boundary layer in more detail in a future paper. 
How much will the boundary layer contribute to the damping (or growth) of the g-modcs? 
Further work should try to explain why the value of Q is so high. Naively, one would expect 
Q ~ (It is fortunate that most of our results are fairly insensitive to the value of a.) 

The nonlinear growth and coupling of the modes must also be investigated. 

This work was supported by grant 075-2002 (Incentives) of Ministerio de Ciencia y 
Tecnologfa, Costa Rica, grant 829-A3-078 of Viccrrcctoria dc Invcstigacion, Univcrsidad de 
Costa Rica, and NASA grant NAS 8-39225 to Gravity Probe B. We are also grateful to the 
Aspen Center for Physics for its support. 
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Table 1. Value of J as a function of j and parity 



j parity J 



even 1/2 

odd 1 

1 even 3/2 
1 odd 2 



-17- 



Table 2. Values of the eigenfrequency a x 10^ for m = g-modes. 



n 




J 




1 


3/2 


2 





2.65 


2.73 


2.75 


1 


2.39 


2.55 


2.61 


2 


2.23 


2.43 


2.51 



Table 3. Values of the modes' fractional radial width Ar/r^ for m = g-modes 



n 




J 




1 


3/2 


2 





0.044 


0.024 


0.019 


1 


0.12 


0.072 


0.055 


2 


0.17 


0.11 


0.083 



Table 4. Harmonic structure of the J —l^n — Q g-modes 



m 











2.65 


1 


0.044 


1 


5.39 


2.03 


0.033 


2 


8.14 


3.07 


0.028 
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Fig. 1. — The horizontal hnes show the radial extent and the frequency of the three lowest 
fundamental {n — 0, J = 1) axial g-modes. 



